function hrun()

  Tfinal = 1;
  Tsteps = 100;
  Nx = 32;
  L = 1;
  Temp = 1;

  params = hparams(Temp, Nx, L, Tfinal, Tsteps);

  params.hydroOrder = 2;

  params.plotHydro = true;

  fname = sprintf('%s', datestr(clock, 'yy-mm-dd HHMM'));



  hinit(params);


  % initial conditions, or replace with your own
  uIn = hic(params);



  % run the ode solver
  [tOut, uOut, udtOut, udtdtOut] = hhydro(params, uIn);

  % save the results
  save(sprintf('%s.mat', fname), 'tOut', 'uOut', 'udtOut', 'udtdtOut', 'params', '-v7.3');






end
